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ABSTRACT 

We calculate the self-consistent response of an axisymmetric galactic disk 
perturbed by an elliptical halo potential of harmonic number m — 2, and obtain 
the net disk ellipticity. Such a potential is commonly expected to arise due to 
galactic tidal encounter and also during the galaxy formation process. The self- 
gravitational potential corresponding to the self-consistent, non- axisymmetric 
density response of the disk is obtained by inversion of Poisson equation for a 
thin disk. This response potential is shown to oppose the perturbation potential, 
because physically the disk self-gravity resists the imposed potential. This results 
in a reduction in the net ellipticity of the perturbation halo potential in the disk 
plane. The reduction factor denoting this decrease is independent of the strength 
of the perturbation potential, and has a typical minimum value of ~ 0.75 — 0.9 
for a wide range of galaxy parameters. The reduction is most important at 
1.4 exponential disk scale lengths and is progressively less so at higher radii. 
For the solar neighborhood region of the Galaxy, the reduction factor is 0.8 . 
Beyond twice the Holmberg radius, the reduction is negligible, and there the 
disk asymmetry in the atomic hydrogen gas traces the true ellipticity of the halo 
potential. The reduction is negligible at all radii for higher harmonics (m > 3) 
of the halo potential. 

On correcting for the negative disk response, the true ellipticity of the halo 
potential for a typical spiral galaxy is shown to be higher by ~ 20% than the 
typical halo ellipticity of < 0.1 deduced in the literature from observations of 
isophotal or kinematical asymmetry of disks. 

Subject headings: galaxies: kinematics and dynamics - - galaxies: spiral - 
galaxies: structure - galaxies:halos - galaxies: ISM 
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1. INTRODUCTION 

It is now realized that the disks of spiral galaxies display a rich variety of non- 
axisymmetry in their light and hence mass distribution. Lopsided galaxies such as M 101 
and NGC 1637 show a global asymmetry, with a much larger spatial extent on one side 
than on the other, as seen optically (Sandage 1961), and particularly strongly in the atomic 
hydrogen gas (Baldwin, Lynden-Bell, & Sancisi 1980). Thus the disk mass distribution in 
these is characterized by m = 1 where m is the harmonic number or the azimuthal number 
of the Fourier component being studied. For the elliptical case (m = 2), similar globally 
asymmetric distributions with a constant phase with radius are not easy to discern. This is 
because it is difficult to separate the projection effects from the intrinsic ellipticity of the 
disk, and hence this needs a careful photometric study of nearly face-on galaxies as done 
for example by Rix & Zaritsky (1995). On the other hand the m = 2 spiral features with 
a phase varying with radius say as in M 81 or M 51 (Sandage 1961) are well studied, and 
indeed are responsible for the 'spiral' nomenclature for these galaxies. These earlier studies 
were done in the blue-band, and hence the features show a much stronger contrast than the 
asymmetry in the underlying stellar mass distribution which is harder to measure. 

Recent near-IR observations allow one to measure the mass asymmetry in the 
underlying, old stellar disk population (Block et al. 1994, Rix & Zaritsky 1995, Zaritsky 
& Rix 1997, Rudnick & Rix 1998, Kornreich, Haynes & Lovelace 1998). Rix & Zaritsky 
(1995) have given a quantitative measure of the amplitudes of the components with various 
harmonic numbers and their radial variations. The amplitudes of the higher harmonics 
(m > 3) are observed to be generally smaller than those for m — 1, 2 or 3. The asymmetry 
at a higher radius beyond the optical disk or the Holmberg radius is better studied using 
the HI as a tracer as shown for the m = 1 case, with mapping (Baldwin et al. 1980) 
and with the global velocity profiles (Richter & Sancisi 1994; Haynes et al. 1998). The 
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disk asymmetry can also be deduced from the kinematic study of the HI velocity field 
observations (Schoenmakers, Franx, & de Zeeuw 1997). 

Despite the small value of the observed disk asymmetry of a few % , it is becoming an 
interesting topic for study since the disk asymmetry is tied into the asymmetry of the halo 
and indeed provides a quantitative measure or diagnostic of the halo asymmetry. This is 
because the global asymmetry in the disk is attributed as a response to the halo distortion 
( e.g., Weinberg 1995, Jog 1997). Thus the halo asymmetry can be deduced from the 
observed isophotal (Rix & Zaritsky 1995), or the kinematic (Franx & de Zeeuw 1992) disk 
asymmetry. Binney (1978) first proposed that a galaxy halo would be non-axisymmetric 
/ triaxial, and studied its effect on the embedded disk. The halo asymmetry could arise 
either due to tidal interactions between galaxies (Wienberg 1995), or the triaxiality of the 
halo could be attributed to the galaxy formation process itself (e.g., Dubinski & Carlberg 
1991; also see Binney 1996, and Rix 1996). Thus the measurement of the disk asymmetry 
allows one to constrain the details of galaxy formation mechanisms. 

Since the disk and halo in a galaxy overlap and interact with each other gravitationally, 
the non-axisymmetry in one structural component will affect that in the other. In the inner 
galaxy, inside of the Holmberg radius, the disk constitutes a significant part of the total 
mass of the galaxy, and the disk self-gravity can result in a substantial decrease in the net 
non-axisymmetry of the halo potential as shown for m — 1 by Jog (1999). Physically this is 
because the disk self-gravity resists change and thus leads to a decrease in the magnitude 
of the net lopsided potential in the galactic plane. In this paper, we study a similar 
self-consistent, negative disk response to higher order halo perturbations (m = 2,3), and 
discuss the implications of our results for observations. A similar opposing disk response for 
m = 2 has been mentioned by Rix (1996), and the decrease in ellipticity of the potential has 
been estimated approximately by Binney (1996). Here we study this effect quantitatively, 
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and self-consistently, and find that the typical reduction in the elliptical potential is ~ 20%. 

The orbits in a m = 2 and m = 3 perturbation potential are calculated using the 
first-order epicyclic theory, and the density response of the disk is obtained (§2). Further, 
the self-gravitational potential corresponding to the self-consistent density response of the 
disk to the halo potential is obtained by applying to m = 2 and 3 the general formalism 
developed by Jog (1999) (§2). The reduction factor for the halo potential due to the 
self-gravity of the disk is obtained for a wide range of galaxy parameters, including for the 
Milky Way, and the 'true' halo asymmetry is obtained (§3). §4 contains a discussion of a 
few general points, and the conclusions from this paper are summarized in §5. 

2. POTENTIAL CORRESPONDING TO DISK RESPONSE 

2.1. Density Response In a Nonaxisymmetric Halo Potential 

We obtain the equations of motion for closed orbits, and the density response of these, 
in an azimuthally symmetric galactic disk perturbed by a non-axisymmetric halo potential, 
and also obtain their relation to the isophotal shapes in an exponential disk (see Appendix 
A). This is analogous to the lopsided case (m = 1) studied by Jog (1997). We use the 
galactic cylindrical co-ordinates (R,<f>,z). 

The unperturbed potential for the axisymmetric galactic disk, ip , at a given radius R 
is chosen to be a logarithmic potential which is applicable for a region of flat rotation, with 
V c being the constant rotational velocity: 

MR) = V c HnR (1) 
The non-axisymmetric perturbation halo potentials ip2 and ipz corresponding to the 
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m = 2 and 3 components are respectively chosen to be: 



ip m {R) = Vc e m cosm <p 



(2) 



where e 2 and €3 are small, constant perturbation parameters, and is the azimuthal angle 
in the plane of the disk. Note that the ellipticity, or the elongation, of the potential in the 
plane is given by 2e 2 and 2e 3 respectively for the m = 2 and 3 cases. The amplitude, V c 2 e m , 
of the perturbation potential (eq. [2]) is assumed to be constant with radius for simplicity, 
as in Franx & de Zeeuw (1992), and Rix & Zaritsky (1995). This assumption is physically 
reasonable for a global distortion of a halo as in a triaxial halo potential (see §1). We 
have also assumed the phase of the perturbation potential to be constant with radius for 
simplicity, as was done by Rix & Zaritsky (1995) for m — 2, and by Jog (1997) for m — 1. 
This is also justified for a global halo distortion assumed. 

The unperturbed surface brightness of a typical galactic disk is observed to have an 
exponential dependence on radius (Freeman 1970). Assuming a constant mass-to-light ratio 
for the disk, this gives the unperturbed mass surface density, /i un , of the stellar disk to be: 



where no is the central extrapolated surface density, and R exp is the scale length of the 
exponential disk. 

For the perturbed case due to the non-axisymmetric potential, the resulting isophote 
will be elongated along the same axis as the perturbed orbit, as argued by Franx & de 
Zeeuw (1992). Since the perturbed orbit, SR, is proportional to cosmcp (see e.g., eq. [A6] 
for m = 2), the resulting effective surface density for the perturbed orbits in an exponential 
disk can be defined (see Rix & Zaritsky 1995, Jog 1997) to be: 



l^un(R) = fioexp[-(R/R exp )} 



(3) 
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fi(R, <p) = /iq exp 
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H exp[ -(R/R exp )] exp[(A m /A ) cosm<f>] 



(4) 



where (e iso ) m is the ellipticity of an isophote at R for m = 2 and m = 3 respectively, and 
A m /A Q is the fractional amplitude of the rath azimuthal Fourier component of the disk 
surface brightness. Note that for small perturbations, this choice reproduces the definition 
of {e iso ) m (eq.[Al2]), as required. The second relation on the r.h.s of equation (4) follows 
from the relation between (ej SO ) m and A m /A (eq. [A13]) for an exponential disk. 

Thus, the change in the surface density, fi reS ponse, resulting from the response of the 
disk to the perturbation halo potential is given by subtracting /x un , the unperturbed surface 
density (eq. [3]), from eq. (4), and for small perturbations we obtain the density response 
for m = 2, 3 to be respectively: 



Thus, the disk response density is linearly proportional to A m /A or to the perturbation 
parameter, e m (see eq. [A15]), as expected from the linear perturbation theory used in 
this paper. Note that the response density is maximum along = 0°, along which the 
magnitude of the perturbation halo potential is also a maximum. This result will be valid 
for any self-gravitating, centrally concentrated realistic disk mass distribution. 
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2.2. Disk Response Potential : \m\ = 2 , 3 Cases 

The self-gravitational potential ip (R, 0, z) for a general, nonaxisymmetric, thin disk 
with a surface density /x (R, , z) was obtained by Jog (1999). This was obtained by 
solving the Poisson equation using the inversion technique involving the Hankel transforms 
of the potential-density pairs. The expression for the nonaxisymmetric potential for the 
thin disk is ( see eq. [22] from Jog 1999): 



°° POO POO 

4>disk (R, 0, z) = —G exp (iriKp) / J m (kR) exp(—k\z\) dk / J m {kR') R' dR' 

, „ Jo Jo 



p2lT 

x / //(#, 0') exp(-im0') dxj>' (6) 
Jo 

where J m (kR) is the cylindrical Bessel function of the first kind, of order m. We apply this 
to the disk response density, l^ r es P onse(R, 0)]m for m = 2 and m = 3 (eq. [5]). The resulting 
potential defines the response potential, (ip r esponse)m = [^response (R, 4>)}m, respectively for 
m = 2 and m — 3. We consider the potential in the plane of the disk, so that z — 0. 

Consider the case \m\ — 2 first. Since the density response is proportional to cos2 ; 
only the values of m = ±2 need to be kept in the integral over on the right hand side of 
equation (6), because only the contribution from these terms is non-zero. Further, since 
cos2 is an even function of 0, and since J 2 (kR) = — J_ 2 (A;i?); hence the terms for m = 2 
and m = — 2 contribute equally to the integral. On substituting from equation (A15) for 
A2/A0 in terms of the expression for the response potential, (ip r esponse)2 , simplifies to: 



r 

(Response) 2 = ~ 2,7rG fX 6 2 COS2 / 

Jo 



J 2 (kR)dk 
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i>oo 

/ J 2 (kR') (1 + R'/R exp ) exp(- 
Jo 



R' 



R, 



)R'dR' 



(7) 



exp 



The second integral over R' can be solved using the relation (6.623.1) from 
Gradshteyn & Reyzhik (1980). The first integral over R' can be simplified by writing 
J 2 (kR') in terms of J\ and J using the standard recursion relation (e.g., Arfken 1970): 
Jo(x) + J2(x) = (2 fx) Ji(x) ; and the resulting two terms can be solved using the relations 
(6.611.1) and (6.623.2) respectively from Gradshteyn & Reyzhik (1980). On substituting 
these in equation (7), and writing x = kR, and setting f(x) = x 2 (R exp /R) 2 + 1, equation 
(7) simplifies to: 



response ) 2 = - 2nG /iQ R exp e 2 cos2 



R 



R, 



exp JO 



J 2 (x)dx 



X 
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_ (Rexp/ ry 

if(x)}V 2 J ~ [/(x)]3/2 



R f c 
- 2nG /2 R exp e 2 cos2 (-|^) 3 J 



J 2 (x)3x 2 dx 

[f\x)Y >/2 



(8) 



Note that (jp r esponse)2 has a sign opposite to the perturbation potential ip 2 (eq. [2]), thus 
the disk response is negative. In the linear regime studied, the magnitude of the response 
potential is proportional to e 2 , the perturbation parameter in ip 2 , and cos2 <p. Similarly, 
the disk response would be negative for any self-gravitating, centrally concentrated disk. 
Following a similar analysis, the response potential ijpresponse) 3, is obtained for m = 3, see 
Appendix B for details. 



Next, define i] m to be a dimensionless quantity as: 
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(Vm) 



response 



) m 



rn 



(9) 



Note that this is independent of the strength of the imposed perturbation potential, and 
depends linearly only on /x R exp /V c 2 . 

The integrals over x in equation (8) when solved analytically give terms involving 
hypergeometric series which need to be calculated numerically. Instead, equation (8) is 
directly solved numerically, and while doing this Jz(x) is obtained using a specific program 
(Press et al. 1986, Chap. 6) that gives a stable value for n > 2. 

Define the dimensionless disk response potential, j m , as: 



In Figure 1, 7 m vs. R/R exp is plotted for a flat rotation curve for m = 2 and 3. First, 
note that, at any radius, the magnitude of 7 m is lower for the higher m value. This follows 
from the form of the response potential (eq. [6]), which involves a double integral over the 
Bessel Function J m {x) which decreases monotonically with m for a particular argument x. 
Second, the maximum of 7 m occurs at a lower radius with increasing m. The maximum 
occurs at a radius = 1.42 R exp for m = 2 and 1.17 R eX p for m — 3. For the corrected 
equations of motion for m — 1 (see Appendix A), the maximum occurs at 1.98 R exp which 
is slightly larger than the radius of 1.4 R exp obtained in Jog (1999). This radial dependence 
is a robust result valid for any exponential disk, and is independent of the actual values of 
/io and R exp . The observational consequences of this result will be discussed in § 4 . 



V c 2 



(10) 



2tvG^oR, 
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3. RESULTS 

3.1. Net Nonaxisymmetric Potential: Self-Consistent Calculation 

A particle in the disk will be affected by the imposed halo potential and also the 
disk response potential. We follow the approach of Jog (1999) which gave a self-consistent 
calculation for m = 1. Thus, for a self-consistent case, the net nonaxisymmetric perturbation 
potential, (ip ne t)m, in the disk plane is given by the sum of the perturbation halo potential 
(eq. [2]), and the self-gravitational potential, (ipresponse)' m i which corresponds to the disk 
response to the net potential, (ip ne t)m- Thus, 

response 

)'m (11) 

In analogy with the disk response to the disk potential alone (eq. [9]), the net self-consistent 
response potential {ipresponse)'m * s gi ven to be: 

response 

Ym = -Vm(^net)m (12) 

On substituting this in equation (11), we get 



(VWt)m = : ~ (13) 
1 +Vm 

Next, define the reduction factor, 5 m , to be: 



*" s TT^ (14) 

Here S(< 1) is the reduction or scaling factor by which the magnitude of ip m is reduced 
due to the self-consistent, negative disk response. Note that since r\ m is a positive definite 
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quantity, hence 8 m < 1. That is, the magnitude of the net perturbation potential is always 
smaller than the magnitude of the imposed halo perturbation potential. For 5 m — 1, there 
is no reduction and the disk response n m = as expected. The reduction factor 5 m at a 
given radius R/R exp is independent of the strength of the perturbation potential and hence 
of e m , and it depends inversely on r\ m and hence inversely on /i R e x P /V c 2 . We will obtain 
the actual values of S m in §3.2 . 5 m will be a minimum at a radius where n m is a maximum, 
that is at 1.42 and 1.17 R/R exp for m = 2 and 3 respectively (see Figure 1). 

3.1.1. Net Asymmetry 

Define the net, self-consistent, perturbation potential, (ip n et)m, in terms of a small 
perturbation parameter (e net ) m to be: 

(^net)m = K 2 (enet)m COS(f) (15) 

Substituting this, and ip m (from eq. [2]) into equation (13), we obtain: 

(e ne t)m = Cm (8)m (16) 

Thus, the parameter e net denoting the strength of the net perturbation potential in the 
galactic disk is reduced compared to e m , the parameter denoting the perturbation halo 
potential, by the reduction factor 5 m . Observations of A m /A , the fractional amplitude of 
the m th azimuthal Fourier component of the surface brightness, will yield the parameter 
{^net)m- Here 2(e net ) m is the net ellipticity of the halo potential. Hence the halo-alone case 
(eq. [A15]) is now modified, and we get the net ellipticity to be: 
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A 2 /A 



W> = T, ' R x (17) 



Similarly, (e net )3 is given by r.h.s. of equation (A7). The values of the true ellipticity 
(2e 2 ) will be obtained in §3.3. 



3.2. Reduction Factor, <5 m 

The value of the reduction factor 5 m (eq. [14]) is obtained numerically and its variation 
with the galaxy morphological type, size, and radius in the galactic disk, and the component 
m, is studied for the classical large or giant spiral galaxies. 



3.2.1. Giant Spiral Galaxies 

The values of the typical disk parameters for the giant spiral galaxies are taken to be: 
Ho, the central surface mass density = 450M Q pc~ 2 ; R exp , the exponential disk scale length 
= 3kpc; and a range of values of the maximum rotation velocity, V c , for a flat rotation 
curve are taken to be equal to 200, 250, and 300 kms' 1 . See Jog (1999) for a discussion 
supporting the choice of these values. Figure 2 contains a plot of the reduction factor, (<5) 2 , 
by which the elliptical halo potential is reduced due to the negative response of the disk (eq. 
[14]), vs. R/R exp . Note that the reduction factor is a minimum at lA2R exp and increases 
thereafter, as expected from equation (14) (see §3.1). The typical minimum value of S lies 
in the range of 0.75 - 0.9 , and 5 is larger for galaxies with a larger value of V c . 

A similar plot of 63 versus radius (not shown here) gives the minimum value of £3 to 
be higher, in the range of 0.83 — 0.93. Thus the reduction due to the disk self-gravity in the 
perturbation halo potential is not important for m — 3, and for the higher values of the 
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harmonic m. Therefore, for m > 3, the observed asymmetry in the disk response represents 
the true halo asymmetry 

For the corrected equations of motion for the lopsided case (m — 1, see Appendix 
A), the resulting minimum in 5 is found to span a range of 0.67 - 0.82 and it occurs at a 
radius of 1.98R exp . These are slightly different from the results of Jog (1999), mainly in 
the peak radius which was earlier obtained to be lAR exp . Note, however, that the correct 
results obtained here are still in good agreement with the observed radial variation in the 
net lopsided distribution, which was a main result of Jog (1999). 

3.2.2. The Milky Way 

Next consider the special case of a giant spiral galaxay, namely, the Milky Way Galaxy. 
We assume a flat rotation curve with V c = 220kms~ 1 , R exp = 3.5/cpc (Binney & Tremaine 
1987), and fi = 450 M G pc~ 2 as discussed above for the giant galaxies. The plots of <5 2 and 
#3 vs. R/R exp are given in Figure 3. Several quantitative resuls follow from this. First, 
the minimum values of 5 2 and 5 3 are 0.75 and 0.88 respectively. Second, 5 2 is 0.79 in the 
solar neighborhood of R = 8.5kpc, which is = 2.43 x R exp . This is almost twice the value 
of 3/7 for the reduction in the ellipticity estimated from an order-of-magnitude calculation 
by Binney (1996). Because of the general formulation in our paper, Figure 2 gives the 
reduction factor as a function of radius for m = 2 for a variety of galaxy parameters. 

Third, if the halo of our Galaxy has an elliptical halo potential, then the disk response 
would reduce this potential at most by a factor of ~ 0.75. Thus for the observed values 
of the disk parameters for the Galaxy and also the other galaxies, while the negative disk 
response cannot be ignored, it can never totally cancel or counteract the imposed elliptical 
halo potential in the disk plane. This is contrary to the suggestion by Binney (1996) that 
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at high enough disk to halo mass ratio, the galaxy could be treated as axisymmetric. 

3.3. Ellipticity of the Halo Potential 

The true ellipticity of the halo potential is an important physical property, possibly 
related to the process of galaxy formation (§1), and attempts have been made in the 
literature to estimate this from the observational data on disks of galaxies. The resulting 
values span a large range. The optical data on the elongation in the disk yield typical 
estimate of ellipticity for the halo of spirals to be 0.1 (Franx & de Zeeuw 1992), while a 
value of 0.045 is obtained by a similar analysis of the near-infrared study of a smaller sample 
of 18 galaxies (Rix & Zaritsky 1995). From detailed kinematical studies, an ellipticity of 
0.1 is obtained for the Milky Way Galaxy (Kuijken & Tremaine 1994), and for NGC 2403 
and NGC 3198 this is estimated respectively to be 0.064 and 0.019 (Schoenmakers et al. 
1997). From the scatter in the Tully-Fisher relation, the halo ellipticity is estimated to be 
0.1 (Franx & de Zeeuw 1992). 

The observations measure the net ellipticity, 2(e net ) 2 , while the true ellipticity of 
the halo potential, given by 2e 2 , is higher by a factor of l/<5 2 (eq.[16]). Since the typical 
minimum value of 5 2 is 0.8 (Figure 2), the true ellipticity is higher by 20 % than the 
measured net value. Thus, for the above observed typical range of net ellipticity of 
0.045 — 0.1 in the literature, the true halo ellipticity is in the range of 0.056 - 0.12 . This 
is an important new physical result from our work. Of course, the above estimates would 
have substantial error bars due to the contamination by the spiral arms or a central bar. 

Note that the reduction factor S 2 — > 1 at large R > 8R exp (Figs. 2-3), or about twice 
the Holmberg radius. Hence the true halo ellipticity can be directly sampled by studying 
the tracer at larger radii namely atomic hydrogen gas. This was done in the plane of IC 
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2006 (Franx, Van Gorkom, & de Zeeuw 1994), and the halo was found to be axisymmetric. 
The ellipticity perpendicular to the plane of the galactic disk could be obtained by studying 
the polar ring galaxies as suggested by Rix (1996). 

4. DISCUSSION 

1. The net nonaxisymmetry in the disk will only manifest beyond the radius where the 
magnitude of the disk repsonse potential is a maximum, and its magnitude will increase 
with radius as seen from the definition of 8 m (§3.1), and this radius is larger for a lower 
m. This indicates the increasing relative dynamical importance of halo over the disk at 
large radii. Therefore, if the halo distortion is constant or increasing with radius, then 
the disk lopsidedness (m — 1) would be seen only in the outer disk while the higher-order 
nonaxisymmetric features could be seen farther in the disk. However, the radial variation 
for m > 2 in the inner/optical region will be affected by spiral arms, and bars. Hence the 
radial variation in m > 2 components cannot be given clearly, in contrast to the m = 1 case 
where a clear minimum radial distance was predicted for the detection of m — 1 global 
features (Jog 1999, also §3.2). 

2. In addition to the disk response to the global perturbation in the halo as studied in 
this paper, there could also be m = 2 or m = 3 modes generated directly in the disk, say 
due to gravitational instabilitites. These would typically have a strong phase variation with 
radius and hence be detected as the standard two-armed or three-armed spiral features 
respectively. Only future detailed simulations of tidal interactions on lines of Wienberg 
(1995) will tell us about the strength as well as the phase and the radial dependence of the 
true halo non-axisymmetry of the various m components. Since there is no Inner Lindblad 
Resonance for m = 1 in a typical galactic a disk (e.g., Block et al. 1994), the m — 1 
component may dominate in the non-linear regime. 
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3. The negative disk response decreases the net asymmetry of the potential in the 
galactic plane, and this would affect the further evolution of the galaxy. This could be one 
reason why numerical simulations with a higher mass concentration in the form of a disk 
show a decrease in the halo ellipticity (Dubinski 1994). Thus the negative feedback due to 
the self-gravity of the disk highlighted in this paper should be included in future studies of 
galaxy evolution. The present paper shows that the disk cannot be treated as a collection 
of massless test particles - doing so would overestimate the disk reponse. 



5. CONCLUSIONS 

We have calculated the self-consistent disk response for an axisymmetric galactic 
disk perturbed by a non-axisymmetric halo potential with elliptical and higher-harmonic 
perturbations (m = 2 and 3) : 

1. The self-gravitational potential of the self-consistent density response of a galactic 
disk is calculated and this is shown to oppose the perturbation potential. Thus, the 
magnitude of the net non-axisymmetric potential in the galactic disk plane is always 
reduced compared to that of the perturbation potential. This reduction is denoted by a 
factor, 5 m , which is found to be independent of the strength of the perturbation potential. 

2. The reduction factor, 62 , is obtained for a wide range of galaxy parameters, including 
for the Milky Way. It has a minimum value of ~ 0.75 — 0.9, which is insensitive to the 
morphological type and size of the galaxy. The reduction is most significant at 1.4 disk 
scale lengths and is less important at higher radii. Beyond twice the Holmberg radius, the 
reduction is negligible and the atomic hydrogen gas can be used to trace the true ellipticity 
of the halo potential. In the solar neighborhood of the Milky Way, the elliptical halo 
potential is decreased by a factor of <5 2 = 0.8 due to the disk self-gravity. The reduction 
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is negligible for the higher harmonics (m > 3) of the halo potential. The asymmetric disk 
response in m > 3 therefore represents the true halo asymmetry. 

3. On correcting for the negative disk response, the true ellipticity of the halo potential 
for a typical spiral galaxy is shown to be higher by ~ 20% than the halo ellipticity 

of ~ 0.5 — 0.1 deduced in the literature from observations of isophotal or kinematical 
asymmetry of disks. 

4. The negative disk response due to the disk self-gravity is shown to be always 
significant in decreasing the strength of the non-axisymmetric halo potential. Hence the 
galactic disk in a realistic galaxy cannot be treated as a collection of massless test particles. 
Yet, the negative disk response in a real galactic disk is found to be never large enough to 
completely obliterate the effect of the halo ellipticity. 

I would like to thank the anonymous referee for comments that led to a clearer 
presentation of material in Section 2, and the Appendix A. 
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A. APPENDIX A - PERTURBED ORBITS AND ISOPHOTAL SHAPES 

We study the orbits in an axisymmetric disk perturbed by a halo potential of harmonic 
numbers m — 2, and 3. We use the cylindrical co-ordinate system (R,(f>), where is the 
azimuthal angle in the galactic plane. Consider a perturbed orbit around the initial circular 
orbit at a radius -Ro, which is given by R = Ro + 5R and = 0o + 5cj). Here 0o = 
where Qo is the circular rotation speed at Ro and is given by the following where ipo is the 
unperturbed potential: 



Ro^o 2 = =k (Al) 

Consider the m = 2 case first. The general perturbation potential is taken to be 
tppert(Ro) cos2<po- Following the procedure for the first-order epicyclic theory as in Jog 
(1997), the general coupled equations of motion for SR and 5(f) are given by equations (4) 
and (5) from Jog (1997) to be: 



-MT=SR [™o + U J - [-^- + k ] co S 20o (A2) 



= — ifc— sm20 ° (A3) 

From the theory of a forced oscillator (e.g., Symon 1960), equation (A2) may be solved 
to yield the following solution for the closed orbits: 



5R = 



2tl>pert(Ro) _|_ dlpp er t I 

Ro dR \ R o 

K 2 - 40, 2 



cos2(po (A4) 



where k is the epicyclic frequency at Rq. 



-20- 



For the present problem, ippert = Vc 2 £2 is the amplitude of the perturbation potential 
■02 defined by equation (2), and the unperturbed disk potential ip is defined by equation 
(1). For these, equation (A4) gives: 

5R = R e 2 cos2(p (A5) 

Thus the net radius is given as: 

R = R + SR = R (1 + e 2 cos 20 o ) (-46) 
Hence Vr, the perturbed velocity along the radial direction, is given as: 

V R =-2V c e 2 sin2(f) (A7) 

On substituting the solution for 5R from equation (A5) in equation (A3) for the 
component of the equation of motion and integrating it we obtain the solution for the 
perturbed azimuthal velocity component R d(5(f))/dt . The net velocity along the azimuthal 
direction is obtained to be: 

V4 = V c + R d(5<f))/dt + tt 5R = V c (1 - 2 e 2 cos 20 o ) (A8) 

Hence the equations of motion for the perturbed, closed orbits in the m = 2 perturbed 
potential are given by equations (A6)-(A8). 

A similar procedure for m = 3 yields the following equations of motion for the 
perturbed, closed orbits: 



- 21 - 



R = Ro[l + (2/7)e 3 co S 30 o ] ; V R = -{6/7)V c e 3 sin3<f> ; = V c [l - (9/7)e 3 co S 30o] (A9) 

The results for the m — 1 (lopsided) case from Jog (1997) are given below for 
comparison. The change in the azimuthal velocity, SV^, was given in Jog (1997) to be equal 
to R d(5(j))/dt. The correct, space-frame velocity (see Schoenmakers 1999), would include 
an additional term = fi SR (see for example the l.h.s. of eq. [A8]). On including this 
correction, the perturbation term in is now changed and has a factor 1 instead of 3 in it, 
while the expressions for R and Vr remain the same as in Jog (1997). The revised equations 
of motion for the m = 1 case are given here for the sake of completeness, and are: 



R = Rq(1 - 2e lC os0 o ) 5 V R = 2V c e lS in<p Q ; = V c (l + e lC os0 o ) (A10) 

Next, we study the resulting isophotal shapes for an exponential disk. This analysis is 
similar to the m = 1 case studied by Jog (1997). For an exponential galactic disk (eq. [3]), 
A m /A , the fractional amplitude of the m th azimuthal Fourier component of the surface 
brightness is obtained to be: 



,./,»/ AR R . . A . 

A m /A = Afi/<fi>=\ - — -—I (All) 

where < \x > is the azimuthal average of the disk density. Here AR/ R is the distortion in 
the isophote and is related to (e iso ) m , the ellipticity of an isophote at R, as follows: 



{e tso ) m EE 1 - (Rmin/Rmax) = 2 (AR/ R) (A12) 

where R m in and R max are the minimum and maximum extents of an isophote respectively. 
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From equations (All) and (A12), we get: 



A m /A = (A13) 

A -tiexp 

We next obtain a relation between the perturbation parameter, e m , and the resulting 
A m /A . Since the orbital velocity changes along the perturbed orbit, the associated surface 
density also changes as a function of the angle . The changes for particles on these orbits 
are governed by the equation of continuity which has the following form in cylindrical 
co-ordinates: 



— [R »(R, 0) V R (</>)] + — [v(R, 0) VM\ = {AU) 

On solving together the equations of the perturbed motion as given by equations 
[A6]-[A8] for m = 2, or equation [A9] for m = 3, or equation [A10] for m — 1; and the 
continuity equation (eq. [A14]); and the equation for the effective surface density (eq. 
[4]); we obtain the relation between (e iso ) m , the ellipticity of an isophote, and e m , the 
perturbation parameter for the potential, at a given radius R. On combining this with 
equation (A13), we get the relation between e m and A m /A , the fractional amplitude of the 
azimuthal Fourier component, for the m — 1, 2, and 3 cases to be respectively: 



Ai/Aq _ ^ = A 2 /A _ ^ = A 3 /A 
(2R/ R exp ) — 1 ' (1 + R/R exp ) 1 + 



This is used to write the disk density response (eq.[5]) and hence the response potential 
^2.2) in terms of e m , the perturbation parameter of the potential. 
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A. APPENDIX B - DISK RESPONSE POTENTIAL : \m\ = 3 CASE 

We obtain the disk response potential, (ipresponse) 3 defined in §2.2 , starting from eq. 
[6]. Following the reasoning as in §2.2 , only the terms for m ± 3 need to be kept, and using 
the expression for A 3 /A (eq.[ A15]), we obtain : 

poo 

response 

) 3 = - 2nG fx e 3 cos3 / J 3 (kR)dk 

Jo 

poo D/ 

x / J 3 (kR')[l + (2/7)(R'/R exp )]exp(——)R'dR' (Bl) 

JO -tiexp 

Next, we use the recursion relations between J 3 and the lower order Bessel functions, J 
and Ji (e.g., Arfken 1970) to simplify the integrals over R', and solve the various terms by 
applications of the relations (6.623.1, .2, and .3) and (6.611.1) from Gradshteyn & Ryzhik 
(1980). Next write x = kR, and define the function f(x) = x 2 (R exp /R) 2 + 1; and obtain: 

poo 

response 

) 3 = -2nGfj J0 Rexpe3Cos3(p (R exp /R) 2 / J 3 (x) dx 

Jo 



-x + 8([/(x)]V2 _ 1} 4(iViW 



(/(X))^ ^(R^/R)^ ^[/(X)] 1 / 2 



—27iGfx R exp e 3 cos30 - / J 3 (x) <ir 



2 

7 Jo 



x 



3x (R exp /R) 2 8(R/R exp ) 2 ^ 1^ 



{f{x)f/ 2 X* V (f(x)Y /2 J X(f(x)) 3/2 



(B2) 
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FIGURE CAPTIONS 

Fig. 1. — The dimensionless self-gravitational potential of the disk response, (i) m versus the 
dimensionless radius R/R exp for the azimuthal wavenumber, m = 2 and 3. The maximum 
occurs at a radius of 1.42 and 1.17 exponential disk radii respectively for m = 2 and 3. 

Fig. 2. — The reduction factor S due to the self-consistent, negative disk response for the 
m = 2 perturbation halo potential versus the radius R/R exp , for giant spiral galaxies, with 
a flat rotation curve with a velocity V c = 200,250, and 300/cms -1 , and R exp = 3 kpc. The 
minimum reduction factor lies in the range of 0.75 - 0.9 , and it always occurs at R/R exp = 
1.42, and 8 increases steadily beyond this radius. 

Fig. 3. — The reduction factor S versus the radius Rj R exp) for the Milky Way, with a flat 
rotation curve with a velocity V c = 220kms~ 1 , and R exp = 3.5 kpc, for m = 2 and 3. The 
minimum reduction factor is 0.75 and 0.88, similar to typical giant spiral galaxies (see Fig. 
2); and occurs at 1.42 and 1.17 R exp for m = 2 and 3 respectively. 



m = 2 Perturbation Potential 




